dc Conductivity of an array of Josephson junctions in the insulating state 
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We consider microscopically low-temperature transport in weakly disordered arrays of Josephson 
junctions in the Coulomb blockade regime. We demonstrate that at sufficiently low temperatures the 
main contribution to the dc conductivity comes from the motion of single-Cooper-pair excitations, 
scattered by irregularities in the array. Being proportional to the concentration of the excitations, the 
conductivity is exponentially small in temperature with the activation energy close to the charging 
energy of a Cooper pair on a superconductive island. Applying a diagrammatic approach to treat the 
disorder potential we calculate the Drude-like conductivity and obtain weak localization corrections. 
At sufBciently low temperatures or strong disorder the Anderson localization of Cooper pairs ensues. 

PACS numbers: 74.81. Fa, 71.30.+h, 73.23.Hk, 74.50. +r 



Artificially fabricated Josephson junction arrays (JJA) 
reveal various fundamental quantum phenomena ranging 
from quantum phase transitions to the motion of single 
charges and vortices (see, e.g. [ij). The flexibility of their 
design and parameters makes JJAs a perfect laboratory 
for study of physics underlying these phenomena. 

Low-temperature transport in the array is determined 
by the ratio of the characteristic Josephson coupling en- 
ergy J to the effective charging energy B /2 of adding one 
Cooper pair on a superconductive island in the array. In 
the two limiting cases, J > Jc and J < Jc, where Jc ~ B, 
the JJA is known to be macroscopically superconduct- 
ing or insulating, respectively, which has been demon- 
strated by a microscopic calculation quite long ago 0. 
Especially interesting is the two-dimensional {2D) case 
when a Josephson- or charging-energy- dominated array 
may undergo respectively a vortex- or charge- unbinding 
Berezinskii-Kosterlitz-Thouless (BKT) [1] transition to a 
normal conducting state. 

Extensive studies of JJA dynamics have been carried 
out since the first array fabrication 43. However, trans- 
port properties of the Coulomb-blockaded JJAs are not 
fully understood yet, although the issue of the conduc- 
tivity of a granulated superconductor in the insulating 
phase has first been addressed already in Ref . [2] . It was 
demonstrated that the frequency-dependent conductiv- 
ity a (w) had sharp peaks at frequencies corresponding 
to the excitation energies of Cooper pairs. However, the 
dc conductivity in the insulating state has not been cal- 
culated explicitly. 

Later, the conductivity of large JJAs has been stud- 
ied close to the superconductor-insulator transition us- 
ing mean-field-type approaches ^] or scaling arguments 
based on the charge- vortex duality @] . In many respects 
the study of transport was phenomenological and, partic- 
ularly, did not account properly for the effects of disorder. 

However, as the relaxation of charge carrying excita- 
tions due to inelastic processes strongly decreases with 
temperature, the disorder must play the major role in 



the low-temperature transport in the insulating phase of 
the array. 

Actually, disorder is intrinsically present in the conven- 
tional tunnelling Hamiltonian describing the coupling be- 
tween superconducting islands, as electrons can hop from 
one island to any state near the Fermi surface in another 
island. Such a disorder determines the conductivity of a 
regular array of normal metal grains (see, e.g., Ref. 0). 
Is the same true also for the insulating state of a regular 
array of superconducting islands or the macroscopic ir- 
regularities in the system should be accounted for? How 
does the conductivity depend on the temperature and on 
the range of electron-electron interactions? 

In the present Letter we study a large JJA deeply in the 
insulating state (J ^ B/2) and address these questions. 
Most of our results are valid for any dimensionality of the 
array, although at some points we restrict ourselves to 
2D. We assume for simplicity that the superconducting 
gap A in a single island is the largest energy scale. 

We calculate the conductivity of the JJA under rather 
general assumptions and demonstrate that the conduc- 
tivity of a regular periodic array without a macroscopic 
disorder remains infinitely large as long as macroscopic 
disorder and inelastic processes are neglected. The dc 
is carried mainly by single-Cooper-pair excitations with 
the charge ±2e. Such bosonic particles move in a regu- 
lar array without being scattered. At the same time, the 
density of the bosons is exponentially low in temperature, 
cx exp(-£'o/T), with Eq close to B/2. 

Macroscopic disorder in the JJA results in the boson 
scattering on irregularities and makes the conductivity 
finite. The conductivity is proportional to the density of 
the bosons, i.e to exp(— i5o/r), in the limit T <^ Ed 



e^T- V(min(T, J)f ex:p{-Ea/T), 



(1) 



where t is the elastic scattering time and Ed is the energy 
of two bosons of opposite charge (boson dipole) located 
on neighboring islands. The energy Ed is either of the 
same order as _Bo or considerably smaller depending on 
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the range of the effective Coulomb interaction between 
the bosons, determined by the capacitive properties of 
the array. Eq. ^ is an analogue of the classical expres- 
sion for the conductivity of free particles. In the limit 
T <^ J the pre-exponential can be evaluated exactly and 
equals Ae'^Tr/'K. 

Remarkably, the description in terms of scattered 
bosons allows one not only to obtain the classical limit, 
Eq. ([T]), but also to describe the quantum interference 
leading to localization effects. In 2D the first weak local- 
ization correction dawL to a takes the form 

SawL ~ -e^T-i min( J, T) exp{-Eo/T) ln{L^/l), (2) 

where / ~ ( Jmin( J, T))2r is the mean free path (mea- 
sured in lattice periods), and is the Cooper-pair de- 
phasing length determined by their recombination, in- 
teraction with phonons, emission of Cooper-pair dipoles, 
etc. In the present paper we assume that the tempera- 
ture is low enough (T < E^) and thus L^:^l. At T < J 
the pre-exponential factor in Eq. ^ equals — 4e^/7r^. At 
low temperatures or strong disorder the correction Sayy l 
becomes comparable with a, which corresponds to the 
strong Anderson localization of the bosons. 

Eqs. H]) and ([2]) are the main results of our paper. The 
conductivity of a JJA in the insulating phase is similar 
to that for electrons in disordered metals but contains in 
addition the activation exponent determining the density 
of the excited Cooper pairs. 

Now we formulate the model and derive the above re- 
sults. We start with a standard effective Hamiltonian 
Ti. describing the motion of Cooper pairs in a JJA 



cos 



0.), (3) 



where the indices i and j label the superconducting is- 
lands, Bij is the inverse capacitance matrix of the array 
in units of (2e)^, (pi and = —id/d(pi are respectively 
the phase of the superconducting order parameter and 
the operator of the number of excess Cooper pairs in the 
island i. Jij is the energy of Josephson coupling between 
neighboring islands i and j, Jij — Jji — J. In these no- 
tations, adding a Cooper pair to the site i requires the 
charging energy Bii/2, which is site-dependent. The av- 
erage value of this energy is 5/2 = {Bii/2). We consider 
the JJA deeply in the insulating phase, J <^ B. 

Irregularities of the array can be described by the fluc- 
tuations 5Bij and 6 Jij. As we show below, current in the 
array is carried by the individual bosons or antibosons de- 
scribing respectively excess Cooper pairs on the islands 
or "Cooper-pair-holes". Random offset charges weakly 
coupled to the array may shift the energy of a Cooper 
pair on an island, and, thus, contribute to the fluctua- 
tions of the coefficients Ba. 

Neglecting the Josephson couplings and disorder one 
obtains a discrete spectrum of the excitation energies of 



the system determined by the eigenvalues of the first term 
in Eq. ([3]). The eigenvalues of the operators fii are in- 
tegers. The ground state corresponds to all rii — 0. All 
excited states are degenerate as long as the system re- 
mains translationally invariant. Of course, no dc current 
can flow through the system in this limit. 

The degeneracy of the excited states is lifted in the 
presence of the Josephson tunnelling Jij between the is- 
lands. The Hamiltonian, Eq. ([3]), is equivalent to a tight- 
binding model for bosons: their states form a band with 
a width proportional to J. As a result, macroscopic dc 
transport is possible. The disorder results in the scatter- 
ing of the bosons inside the bands and leads to the finite 
conductivity, Eqs. H]) and ([2|). 

At low temperatures the dc current is carried by bosons 
with the charge 2e and antibosons with the charge — 2e. 
The conductivity is dominated by the lowest energy 
bands of bosonic and antibosonic states. These two 
bands are located near the energy B/2 8^. 

In order to calculate the conductivity at low tempera- 
tures we may thus consider a reduced Hilbert space: the 
i-th island has only three quantum states: |0)i, zero ex- 
cess Cooper pairs on it, and | ± 1)^, one Cooper pair (an- 
tipair) . It is convenient to rewrite the Hamiltonian ([3|) in 
this space in terms of pseudospin operators , and 
Sf: St\0), = V2|±l)„ Sf\Tl)^ = V2|0)„ Sf\±l)^ = 0, 
Si\0)i — 0, Sf\±l)i — ±|±l)i corresponding to the pseu- 
dospin Si ~ 1. The reduced Hamiltonian 

fire, = i ^ B,,SIS] -\Y^ J^^SfSj (4) 

i,0 



2^^ 



is equivalent to an anisotropic Heisenberg spin-1 model. 
The pseudospin operators obey the conventional commu- 
tation relations 



= 25ijSl 



Sa . s ^ 



(5) 



We assume for simplicity that the islands in the array 
form a square lattice and begin with calculating the ex- 
citation spectrum in an ideal JJA without disorder. The 
states of the Hamiltonian TLred can be classified by the 
5^-projection of the total spin. The ground state corre- 
sponds to — for all i. 

In order to calculate the conductivity we consider 
states corresponding to a single boson or antiboson in 
the array {S^ = ±1). For these states the eigenenergy of 
the first term of the Hamiltonian Tired, Eq. (|4]), equals 
_B/2. In the limit J <^ B, we approximate the eigenfunc- 
tion of Tired for 5*^ = ±1 by a plane wave 



N- 



(6) 



where N is the number of the islands in the array. The 
corresponding excitation spectrum takes the form 



£;(k) ^ B/2~ 2Jcosk^ - 2 J cos ky 



(7) 
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where k = {kx, ky). Thus, the excitation spectrum has a 
narrow band of the width 8 J separated from the ground 
state by the gap 



En = B/2 - 4J. 



(8) 



As long as the gap significantly exceeds the temper- 
ature the density of the bosons is exponentially small. 
In this limit the interaction between them can be ne- 
glected and we can describe the system in terms of a 
single-particle tight-binding Hamiltonian with the spec- 
trum given by Eq. ([7]). Disoder manifests itself in the 
model through the fluctuating parts SB and 5 J of the 
parameters B and J. 

The conductance of the array can be calculated using 
the standard Kubo linear-response theory. The operator 
lij of the current between the z-th and j-th islands reads 



sts: 



(9) 



Its expectation value /y (w) can be expressed through 
the retarded correlation function of currents Iiij^ki{^^)'- 



n 



ij,kl 



(/,,(r)/fc,(0))e'""^dT 



, (10) 



(ki) 



ny-fci(o))(EifeO, (11) 



where E is the electric field, 1^;- vector connecting islands 
k and I. 

Depending on whether the length of the array is 
smaller or larger than the mean free path I, the transport 
in the sample is respectively ballistic or diffusive. 

Using Eqs. (fTU]) and pTjl in the ballistic limit we find 
the conductance of an x Ly rectangular array 

G = Ly sinh (2 J/T) Iq (2 J/T) exp (-Eq/T) . (12) 

vr 

Here Iq is the modified Bessel function. Of course, at I ^ 
Lx Eq. matches the diffusive conductance aLy/Lx 
with tr given by Eq. ([1]). 

Let us note now that in the low-temperature limit the 
array resembles a conventional semiconductor. Indeed, 
at low density of the bosons the particle statistics is not 
important. The bosons thus can be considered as doubly 
charged electrons thermally activated to the conduction 
band of the width 8J, the latter being separated from 
the valence band of the semiconductor by the gap Eq, 
Eq. ([8]). Hence, one can evaluate the conductivity in the 
diffusive regime using the standard diagrammatic tech- 
nique _l9|- We assume that the disorder is weak enough 
for the elastic scattering time r to exceed the inverse 
characteristic kinetic energy min( J, T) of the bosons in 
the conduction band: 



To average over the disorder we assume for simplicity 
a Gaussian distribution for deviations SBa and 6Jij from 
the average values B and J with correlations 



(SB.iSBjj) = fiStj, (SJijSJki) = f2iS'i.k5ji + SaSjk). 

Our results Eqs. ([T]) and ^ with appropriate t apply 
nevertheless for arbitrary not short-correlated distribu- 
tions of the fluctuations. 

The basic element of the perturbation theory is the 
contraction rule for the effective disorder potential V 

(VkipiVk.p,) = (27r)27(ki+k2,ki-p2) 

(5(ki +k2 -pi -P2), (14) 
7(k+,k_) = /l/4 + /2 5](£ 



where 1^ (i = 1 . . . 4) are the unit vectors connecting an 
island with its nearest neighbors. 

The standard procedure (see e.g. Ref. [9]) of evaluation 
of r in the limit /c <C 1 under the condition (fT3|) gives 



(/i/8 + 4/2)/J. 



(15) 



The scattering time t remains of the same order of mag- 
nitude at arbitrary momentum fc ~ 1. 
We use the Kubo-Greenwood formula 



2(2e)2 f dp f de , , , 
(z;„G^(p,e)z;0G«(p,£)), (16) 



r(min(J,r)) > 1. 



(13) 



where v = 2J(sinA;a;,sinfcj,) is the velocity of the boson, 
n{e) ~ exp(— e/T) is the Boltzman distribution of the 
excitations in the conduction band and G^'^{k.,e) are 
retarded (advanced) Green's functions of non-interacting 
bosons, (G^-^(k, e)) = (e - £;(k) ±i/{2T))-\ Using the 
condition we come to Eq. 

We emphasize that Eq. ^ has been obtained in the 
standard scheme neglecting diagrams with crossing im- 
purity lines [9^1. In disordered metals this approach is 
applicable in the limit epT 3> 1, where ep is the Fermi 
energy. Here the role of large parameter is played by 
T(min(J,T)) [cf. Eq. 1^]. 

Next we calculate the weak localization correction 
SawL to the conductivity. Again, the condition (fT5|) 
allows us to repeat the summation of the diagrams of 
Ref. [l3l and to arrive at Eq. 

The limiting case C/Cq ^ 1, where C is the mutual 
capacitance of neighboring islands and Cq is the self- 
capacitance, often corresponds to the experimental sit- 
uation (see, e.g., Ref. (I4I) and is especially interesting 
from the theoretical point of view. In this case the charg- 
ing energy Eq of a single boson significantly exceeds the 
dipole energy Ed, Eo/Ed ~ ln(C/Co) > 1. 

Such bosons resemble vortices in superconductors [l3| . 
A finite ratio C/Cq determines the scale of the interaction 
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of the bosons and plays the same role as the penetration 
depth cutting the logarithmic interaction of vortices in 
superconductors. In the limit C/Cq oo the energy 
of the bosons logarithmically diverges with the size L of 
the sample. Using the analogy with the vortices one can 
expect in this limit the BKT transition, Ref. 0], with the 
critical temperature Tk of order of the dipole energy Ed- 

Properties of the system of the bosons are similar to 
those of the system of vortices in conventional 2D su- 
perconductors. In the BKT scenario single vortices do 
not enter the system below Tr- and the system is a su- 
perfluid. In the JJA considered here, there are no single 
bosons in the limit C/Cq — s- oo and the conductivity van- 
ishes. However, the BKT transition in superconductors 
is known to smear because the energy of a single vortex 
is finite due to the finite penetration depth, which makes 
vortices itinerant resulting in a finite resistivity. In the 
JJA the energy of the bosons is finite due to the finite 
C /Co and their motion makes the electric transport pos- 
sible. As the density of the bosons is proportional to 
exp(— i?o/T), so is the conductivity, Eq. ([T]). At temper- 
atures r <C -Ed ~ Tk one can neglect the presence of 
dipole excitations in the array which could lead to an ad- 
ditional scattering or relaxation of the single bosons and 
to a screening of Coulomb interaction. 

At T > Tk single vortices in 2D superconductors exist 
due to the entropy contribution, Ref. Q, and the super- 
fluidity is destroyed. 2D superconductors, strictly speak- 
ing, do not exist. Thus, as the superfluidity below Tk is 
either absent, one obtains a crossover from the finite ex- 
ponentially small resistivity contributed by the motion of 
single vortices to the resistivity of normal metals. Anal- 
ogously, in a JJA one can expect a crossover from the 
exponentially low conductivity at T Ed, Eq. ([T]), to 
a temperature- independent conductivity at T > Ed- In 
other words, the Coulomb blockade effects are impor- 
tant at low temperature T < Ed but can be neglected 
dX T > Ed- Interference of the bosons results in the ex- 
istence of one more temperature region where effects of 
localization can play an important role and the Coulomb 
blockade is further enforced by the Anderson localization. 

Recently, a similar model of a 2D JJA (but without 
disorder) has been suggested in Refs. 14] and ,J^] to 
describe the ex per iments on strongly disordered super- 



conductors 15|, ll6| . The authors obtained an exponen- 



tial behavior of the conductivity with the activation gap 
B/2 for T > Ed and a double exponential behavior for 
T <^ Ed- Clearly, our results, Eqs. (P), ©, and the ab- 
sence of the Coulomb blockade at T > Ed are in a strong 
disagreement with those findings. A detailed criticism 
of the theoretical treatment of Refs. 14 1 and [l3| can 



In conclusion, we calculated the dc conductivity of 
a large Coulomb-blockaded Josephson junction array at 
low temperatures and demonstrated that it is determined 
by the thermally activated single-Cooper-pair excitations 
on the superconducting islands. In the absence of macro- 
scopic structural disorder in the array the transport is 
ballistic. In the presence of sufficiently weak disorder the 
conductivity in a certain temperature range is described 
by the Drude-type formula multiplied by the activation 
exponent exp(— _Eo/T'), where Eq is the Coulomb energy 
of the single-Cooper-pair excitation. Lowering the tem- 
perature below some characteristic value results in even 
faster decrease of the conductivity due to the Anderson 
localization of Cooper pairs. 

The authors thank M. Yu. Kharitonov for discussions. 
The work has been financially supported by SFB Tran- 
sregio 12, SFB 491, and US DOE contract No. DE-AC02- 
06CH11357. 



be found in Ref. [l3|. A comparison of our results with 
the results of relevant experimental works will be pre- 
sented elsewhere. Also, it would be interesting to clarify 
the question of the many-body localization 11] of the 
bosons in the system under consideration. 
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